Complete asymptotic expansions for eigenvalues of 
Dirichlet Laplacian in thin three-dimensional rods 

D. Borisov 1 , *and G. Cardone 2 



1 ) Bashkir State Pedagogical University, October Revolution 
St. 3a, 450000 Ufa, Russia 
F ^ . 2) University of Sannio, Department of Engineering, 

P-h ! Corso Garibaldi, 107, 82100 Benevento, Italy 

borisovdiOyandex . ru, giuseppe . cardone@unisannio . it 

Abstract 

I We consider Dirichlet Laplacian in a thin curved three-dimensional rod. 

The rod is finite. Its cross-section is constant and small, and rotates along 
the reference curve in an arbitrary way. We find a two-parametric set of 
the eigenvalues of such operator and construct their complete asymptotic 
expansions. We show that this two-parametric set contains any prescribed 
number of the first eigenvalues of the considered operator. We obtain the 
complete asymptotic expansions for the eigenfunctions associated with these 
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■ first eigenvalues 

o ' 
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O ■ Introduction 

> '■ 

The asymptotics of the spectra of elliptic operators in thin domains were studied by 
many authors, see, for instance, [2]- [ID], [12], [H], and the references therein. There 
are two types of thin domains usually considered, namely, thin rods and thin plates. 
Both types were considered in the book |12j . The eigenvalues for elliptic operators 
with Neumann boundary condition on the lateral surface of the rods and on the 
bases of the plates were studied. The asymptotic expansions for the eigenvalues 
and the eigenfunctions were constructed and justified. We also mention the survey 
[D] on thin rods. 

A thin two-dimensional domain formed by two different thin rectangles was 
studied in [H], i.e., the boundary of the domain was non-smooth. The operator 
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considered was Laplacian subject to Neumann condition on the bases and to Dirich- 
let one on the lateral boundary. The main results were the asymptotic expansions 
for the eigenvalues remaining bounded as the width of the domain tends to zero, 
and the asymptotics for the associated eigenfunctions. These asymptotic expan- 
sions were rigorously justified. Laplacian in a thin two-dimensional domain was 
also considered in [TJ. The domain had a variable width with the unique point 
of maximum. The uniform resolvent convergence was established and two-terms 
asymptotics for the eigenvalues were obtained, as well as convergence theorems for 
the associated eigenfunctions. In [8] these results were extended for an infinite thin 
strip under similar conditions for the width. We also mention the paper [10] , where 
a thin strip (finite or infinite) was considered with Neumann condition on the up- 
per boundary and with Dirichlet condition on the lower boundary. Here two-terms 
asymptotics for the first eigenvalues were constructed. The case of a curved infinite 
strip was also studied in [5] , where the number of the discrete eigenvalues below the 
essential spectrum was estimated. The results of [7] were also extended in [2]. Here 
a two-parametric set of the eigenvalues was found and their complete asymptotic 
expansions were constructed. 

A finite three-dimensional rod was considered in [3J. The cross section was 
supposed to be constant and to rotate along the reference curve in an arbitrary way. 
Two-terms asymptotics for the first eigenvalues were constructed and convergence 
theorems for the associated eigenfunctions were established. Similar results were 
obtained in [6] for a tube in a space of arbitrary dimension. An infinite three- 
dimensional thin tube with a round cross section was studied in [5]. The number 
of the discrete eigenvalues below the essential spectrum was estimated and their 
complete asymptotic expansions were constructed. We also mention the paper [I], 
where a multi-dimensional thin cylinder with distorted ends was considered. The 
operator studied was Laplacian in such domain subject to Dirichlet condition on the 
lateral surface and to Neumann one on the distorted ends. The attention was paid to 
the localization effect of some eigenfunctions at the distorted ends. The asymptotics 
for these eigenfunctions and the corresponding eigenvalues were constructed. 

In this paper we extend the results of [3]. We again consider Dirichlet Lapla- 
cian in a curved thin rod. The cross section of the rod is a fixed domain, which 
rotates along the reference curve in an arbitrary way. In what follows this opera- 
tor, its eigenvalues and eigenfunctions are referred to as the perturbed ones. We 
find a two-parametric set of the perturbed eigenvalues and construct their com- 
plete asymptotic expansions. The eigenvalues are indexed by the first two terms 
of their asymptotic expansions. Namely, the leading terms are determined by the 
eigenvalues of Dirichlet Laplacian on the cross-section of the rod. Each of the lead- 
ing terms determines a certain operator on the reference curve, and its eigenvalues 
are the next-to-leading terms of the aforementioned asymptotic expansions for the 
perturbed eigenvalues. It is convenient to group the perturbed eigenvalues into a 
countable set of series, where each series consists of the perturbed eigenvalues with 
the same leading term in the asymptotic expansions. We show that the series as- 
sociated with the smallest leading term contains any prescribed number of the first 
eigenvalues of the perturbed operator provided the rod is thin enough. We prove 
that these eigenvalues are simple and construct complete asymptotic expansions for 
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the associated eigenf unctions. 

In conclusion to this section, we describe briefly the contents of the paper. The 
next section contains the description of the problem and the main results. In the 
third section we introduce a change of variables required for the constructing the 
asymptotic expansions. In the fourth section we select the aforementioned two- 
parametric series of the eigenvalues and construct their asymptotic expansions. In 
the last fifth section we describe the first eigenvalues of the perturbed operator and 
give the asymptotic expansions for the associated eigenfunctions. 

1 Formulation of the problem and the main re- 
sults 

Let x = (xi,X2,xs) be Cartesian coordinates in R 3 , 7 be a finite infinitely dif- 
ferentiate curve in R 3 without self-intersections. By s and sq we denote the arc 
length and the length of 7, s G [0, so]. We parameterize 7 by its arc length, and 
r = r(s) is the infinitely differentiable vector describing 7. The tangential vector 
of 7 is indicated by r = t(s). By 77 = 77(5) we denote an infinitely differentiable in 
s G [0, s ] unit vector defined on 7 being orthogonal to r(s) for all s G [0, So]. We 
let (3(s) := [t(s),t](s)}. It is clear that (3(s) is infinitely differentiable in s G [0,s ], 
and (r, 77, (3) is an orthonormalized frame on 7. One of the possible choices of 77 is 

r](s) := cos a(s)n(s) + sin a(s)b(s), (1.1) 

where n = n(s) and b = b(s) are the normal and binormal vectors of 7, and 
a(s) G C°°[0, So] is an arbitrary function describing how our frame rotates with 
respect to the Frenet one. It follows from (II. ip that 

(3(s) := — sin a(s)n(s) + cos a(s)b(s). (1.2) 

Although this formula and (11. ip could be an appropriate definition of 77 and (3, 
we do not use this way. The reason is that the Frenet frame does not exists for 
all smooth curves, since the normal vector can be undefined at the points where 
r"(s) = 0. 

By u we indicate a bounded domain in R 2 with infinitely smooth boundary, and 
the symbol e stands for a small positive parameter. We introduce a thin curved rod 

as 

n £ := {x G R 3 : x = r(s) + e£ 2 r}(s) + e&/3(s), s G (0, s ), (6, £3) e 

Since the curve 7 is smooth and not self-intersecting, the rod Q e has no self- 
intersections for e small enough. Hereinafter the parameter e is assumed to be 
chosen exactly in such way. 

The main object of our study is the spectrum of the Dirichlet Laplacian in 
L 2 (fl e ), and this operator is denoted by "H £ . We introduce this operator rigourously 
as the Friedrichs extension of — on C£°(Q e ). For each e > the operator % E has 
a compact resolvent and its spectrum is thus purely discrete. The aim of this paper 
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is to construct the complete asymptotic expansions for the eigenvalues of % E . We 
also observe that the eigenvalues of % E can be equivalently regarded as those of the 
boundary value problem 

-A^ e = A e ^ e in Q e , iJ £ = on dQ e , ij} e 6 Wl(VL e ). 

In order to formulate the main results we need to introduce additional notations. 
Let Wf (u;) be the subspace of W^ioj) consisting of the functions vanishing on dui. 
In the same way we introduce the space W^ 2 2 (0, so). By S we indicate Dirichlet 
Laplacian in L 2 ((jj) with W^ofw) as the domain. This operator is self-adjoint. Let 
A n be the eigenvalues of S arranged in the ascending order with the multiplicities 
taken into account, 

Ai < A 2 < A 3 < . . . < A n < . . . 

By (f) n we denote the associated eigenf unctions orthonormalized in L 2 {oS). By the 
smoothness improving theorems [TU Ch. IV, Sec. 2.3] the functions <j) n are infinitely 
differentiable in U. 

Let an eigenvalue A n be simple. Denote 

n:=^—-^ 2 — , C n {u):= / \n<j) n \ 2 di. (1.3) 
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It is straightforward to check that 

t' = KxTJ - T}' — -K\T + K 3 (3, (3' = K 2 T - K 3 T], (1.4) 

where (s) G C°°[0,So] are certain functions characterizing the geometric 

properties of 7 and rotation of rj. By C n we denote the operator 

in L 2 (0, s ) with the domain W / 2 2 (0,s ). The operator C n is self-adjoint. Since 
the operator C n is one-dimensional, by Cauchy theorem we conclude that all its 
eigenvalues are simple. We indicate these eigenvalues by A n ' m ' ) , m = 1,2, . . . and 
arrange them in the ascending order with the multiplicities taken into account. Let 
^(n,m) ^ associated eigenfunctions orthonormalized in L 2 (0, s ). 
If one chooses r/ and (3 in accordance with (11.11) . fll.2p . it yields 

K\ — k cos a, k 2 = k sin a, k 3 = a' + x, 

where k = k(s) and x = x(s) are the curvature and torsion of 7, and k, k G 
C°°[0, s ]. In this case the operator C n becomes 

d " +C n (a;)(a'(s) + x(s)) 2 k2(s) 



Our first result gives the complete asymptotic expansions for the eigenvalues of 
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Theorem 1.1. Let X n be a simple eigenvalue of S. Then there exists a two- 
parametric set of the eigenvalues Xe ,m ' ofH e with the asymptotics 

oo 

X (n,m) = £ -2 Xn + X (n,m) + ^gi^K"^ ^ 

i=l 

where 

Q M ■= (2Ai"' m) - (2C n (u) - i) 4) q + i0 + I^TZq, (1.6) 

q = q( S ,0 := - K 2 (s)^ = := 

The remaining coefficients in the asymptotic series U.5\) are given by the formulas 
in Lemma\JJl 

The series (11.51) are the asymptotic ones and we know nothing on their conver- 
gence. 

We observe that the set of the eigenvalues described in Theorem 11.11 is two- 
parametric and is indexed by n and m. Given n, the eigenvalues \^ n ' mS> (e) form a 
series with the same leading term. We stress that Theorem 11.11 does not imply that 
the eigenvalues Ae"' m ^ form the whole set of the eigenvalues of 1-L e . The first reason 
is the assumption on the simplicity of A„. And even without this assumption it 
is an independent problem to find out whether the eigenvalues \^ n ' m ^ are the only 
possible ones or not. 

We make the assumption that X n is simple in order to simplify the calculations 
in the formal constructing of the asymptotic expansion, see Sec. 3. If A n is multiply, 
it is also possible to construct the asymptotic expansions, but the formal construct- 
ing becomes more complicated and requires some additional careful calculations. 
Another interesting issue is the multiplicities of the perturbed eigenvalues corre- 
sponding to a multiple eigenvalue A n . In view of these issues we regard the case of 
multiple eigenvalue A n as an independent problem, which we postpone for another 
article. 

One more interesting question is on the asymptotic expansions for the eigen- 
functions associated with \i n ' m \ As usually, to justify such asymptotic expansions, 
one has to know lower bounds for the distances between the perturbed eigenvalues. 
The structure of the eigenvalues Ai"' m ^ is such that it is rather difficult to obtain 
such bounds once the eigenvalues are bigger than e~ 2 A2- If we consider only the 
first eigenvalues of % e lying between £" 2 Ai and £~ 2 A 2 , it is possible to prove the 
mentioned lower bounds and to obtain the asymptotic expansions for the associated 
eigenfunctions. This is our second main result. Before formulating it, we introduce 
two additional notations, 

Q := {(s,f) : < s < s ,£ e w}, fl (1) :={(s,():i<s<So-^ew}, tG(0,s / 2 )- 
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Theorem 1.2. Given any M ^ 1, i/iere exists e = £o(M) > swc/i that for 
all e < Eq the first M eigenvalues ofH e are X^e ,m \ m — 1, . . . , M, satisfying U.5\) . 
These eigenvalues are simple and the associated eigenfunctions have the asymptotics 

oo 
i=l 

where the coefficients of the series are given in Lemma \3.1[ The asymptotics hold 
true in W^ity-norm and C k (Q^)-norms for all k ^ ; t e (0,So/2). 

The results of [3] consist of the two-term asymptotics for A^ 1 and the lead- 
ing term in the asymptotics for the eigenfunctions ipe 1 . The asymptotics for the 
eigenfunctions were obtained in /^(fi). Theorem 11.21 extends these results in two 
directions. First, it gives the complete asymptotic expansions. Second, the asymp- 
totics of the eigenfunctions are given in a stronger norm. One more extension is 
provided by Theorem 11.11 Namely, in addition to the first series of the eigenvalues 
Xi 1 ' 171 ^ described in [3], we describe a countable set of similar series 

One more difference from [3] is the technique employed. The study in [3] was 
based on T-convergence of certain functionals. Our approach consists of two main 
steps. The first step is the formal constructing of the asymptotic expansions for the 
eigenvalues and the eigenfunctions by the multiscale method [T]. The second step 
is an estimating of error terms by a result from spectral perturbation theory, see 
Lms 12,13]. 

2 Change of variables 

In this section we transform the operator % e to another one which will be more 
convenient in proving Theorems 11.11 H-21 

Let y = (2/2,2/3) be Cartesian coordinates in the plane spanned over 77 and j3 
with the axes along these vectors and so that the variable y<i corresponds to 77 and 
the variable j/3 does to (3. We first pass to the variables (s, y) and the domain Q e is 
mapped onto 

tt £ := {(s,y) : < s < So,e^y E u}. 

At the second step we rescale the variables y passing to variables (s, £), where 
£ = (£2, £3) = £ -1 2/- Then the domain fl e is mapped onto Q. 

We introduce the operator describing the passing to the variables (s, y) as 

(Uu)(s,y) = u(x(s,y)), {U^u)^) = u(s(x),y(x)). 

We let 

H £ := -pUHjA' 1 , p = p0> v) '= 1 - q(s, y), qO, y) = ^i{s)y 2 - ^O)^. 

If X e and ip e are an eigenvalue and an associated eigenfunction of 7i e , it is clear that 
the function Uip £ solves the equations 

UU £ W x U^ e = X e Ui> e , H £ We = KvW e . (2.1) 
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Let us obtain the differential expression for % E . Taking into account (jl-4j) . and 
differentiating the identity 

x = r(s) + y 2 r]{s) + y 3 f3{s) 

with respect to s, y 2 , and y 3 , we obtain 

dx dx dx 

— =pt(s) - KsysVis) + n 3 y 2 (3(s), — = 77(5), — = (3(s). 

os dy 2 dy 3 

Thus, the derivatives with respect to x and (s, y) are related by the identity 

V (Si2/) = PV X , 
where P is the matrix with the rows 

P:= 77 

It is easy to check that 

detP = p, = P _1 V( S; j,), P^ 1 = (p^V k 3 2/3P _1 t + T7 - « 3 ?/2P _1 t + /3) , 

(2.2) 

where the vectors in the definition of P 1 are treated as columns. Since e 'i/Gw, 
we have y = 0{e) and hence q(s,y) = 0{e). It yields that the function p(s,y) is 
strictly positive for the considered values of y. 
By (I2T2]) for each ui,u 2 G C£°(f2 £ ) we have 

(^wi,«2) ia (n e ) = (P w ^ w ~ ln i' u 2) L2( ^) = iH £ U~ l u 1 ,U' 1 u 2 )L 2 {n E ) 

= (V x W _1 mi, V x W _1 M 2 )L 2 (n £ ) = (P _1 V( a) j / )Wi,pP -1 V( ai j / )ii 2 )x, 2 (n e ) 
= -(div (S)2/) p (p-^P^V^Ui, U2) ia( n e) . 

and therefore 

"H £ := -div( S)J/ ) AV( s , y ) (2.3) 

/ p _1 /^p" 1 -K32/2P _1 
A = ( A ij)i,j=ij = K-mp" 1 p + ^Ip -1 -^s^ysp -1 
\-K3y2P -1 -«iz/22/3P _1 p + ^Ip -1 , 

Now we pass to the variables £. It leads us to a finally transformed operator 

X) f d ( E ) d 9 ( E ) d _j A 9 ( E ) 9 _ 2 A d .(e) 9 \ 

(2.4) 

/ P7 1 ^a&Pj 1 -^s^Pj 1 \ 

A (£) := (4% =13 = [sk&v- 1 p £ + £ 2 Klefp £ - 1 -eWaP," 1 , (2-5) 

\-e«36p e 1 -£ 2 «l6C3P e 1 P e + £ 2 ^lHVe V 

where p e (s, £) := 1 — eq(s,£), and the operator "H £ is considered in L 2 (Q) as the 
Friedrichs extension from C^°(Q). The eigenvalue equation (12.11) casts into the form 

U e ^ e = A £ p £ ^ £ , (2.6) 

where we redenoted (Z/0 e )(s, e£) by ip £ . 
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3 Proof of Theorem 11.1 



In this section we prove Theorem 11.11 The proof is divided into two parts, the first 
being devoted to the formal constructing of the asymptotic expansions. The second 
part consists in proving the existence of the mentioned two-parametric set of the 
eigenvalues and in the justification of the formal asymptotic expansions for these 
eigenvalues, i.e., establishing estimates for the error terms. 

We construct the asymptotic expansions for the eigenvalues and the associated 
eigenfunctions as the series 

oo 

A (n,m) = e i X < [ l ' m \ (3.1) 

i=-2 

oo 

4^\x) = J2e'4 n ' m \s,0. (3.2) 

i=0 

The aim of the formal constructing is to determine the coefficients of these series. 

ij 



We expand the functions A-f in the powers of e, 



4 ) = E^ J) ' (3.3) 





k=0 








4 n) 


= q fc , 






k ^ 0, 


4 12) 


= 0, 


A k 12) 




fc ^ 1, 


.(13) 
^0 


= 0, 


4 13) 


= -«36q fc_1 , 


fc ^ 1, 


.(22) 
^0 


= 1, 




_ n 4(22) _ 2t2 fc-2 


fc ^ 2, 


.(23) 
^0 


= 4 2 3) = o, 


/l(23) _ 2t t n k-2 


fc ^ 2, 


.(33) 
^0 


= i, 


Af 3) 


_ n 4(33) _ 2a2 k-2 
— ^ n k — K 3S2q 5 


k ^ 2. 



(3.4) 



Hereinafter q = q(s,£), if else is not specified. 

We substitute ([23]), lEOl . CT) . lEOl . (ET3j) into (J232) and equate the coefficients 
at the same powers of e. Calculating the coefficient at e l ~ 2 , i ^ 0, we obtain 



(3.5) 



(_a € - aLt } )# ,m) = &4-f + F i n,m) in n > 

^ m ) = o on on, 

F (n,m) ;= _ A («.™) q) ^(«J») > (3. 6) 

d d d d 

1=2 Sl Z=2 S U=2 SI s ' 



=|-q^ 2 |- + *W" 2 |- + |- W"^ + k 2 ^- 2 ^, j £ 2. 
as as os os 
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Consider the problem (13. 5 j) for z = 0. It is clear that its solution can be chosen 

as 

4' Vn \*,t) = *fr n \*)Mt), ^-2 m) =^ (3.8) 
where the function xfr^ 1 '" 1 ^ is unknown and should satisfy the boundary conditions 

^ m) (0) = ^ m) (s ) = 0. (3.9) 
Taking into account (13.81) and an obvious formula 

(9 9 \ 

{T X - A n q)0 re = - ( K!— - K 2 gT- J (f> n , (3.10) 

we write the problem (13.51) for i — 1, 

(-A ? -A n )^i '=-%' M KX—-K2-Q- \ 0n + Ali V„ m fi, 

vi n ' m) =o on an. 

(3.11) 

Employing the eigenvalue equation for <f) n , by direct calculations we check that 

(-A € - \ n )q<f> n = -2 f«i J- - « 2 J-J n . (3.12) 

Hence, the problem (13. lip is solvable for 

X^i m) = (3.13) 

with a solution given by the identity 

#' m) (*,0 = \^ m) (s)MZMs,0 + *i"' m) ( S )0 n (O, (3.14) 
where the function \|/^ n ' m - ) is unknown and should satisfy the boundary conditions 

The above obtained formula for \^ m ^ and that for Aq™'"^ in (I3.8P prove the formulas 
for the first terms in (II. 5ft . 

We substitute the formulas 1EO)]) . fl3~7D . fl3T8|) . fl3~T0l) . fl3~T3l) into the problem 
(13. 5p for i = 2 to obtain 

(-A £ -A n )4 n ' m) = A^ m VS n ' m) + ^ 2 (n ' m) m 0, 

(3.15) 



# im) = o on on, 



9 9 



Ft m) = '(Fl ~ Anq)q^ n ' m) + ~ ^' m) U^- " «2^-J 0n- (3.16) 

Taking into account the eigenvalue equation for <fr n , by direct calculations we check 

(9 9 \ 

(J"i - A n q)q0 n = -3q ( ~ K2 ~9i) ^ n ~ ^ + K ^ n ' ( 3 - 17 ) 

9 



Hence, we can rewrite the formula for F^ 1 '™^ as follows, 

T-,(n,m) 9- / (n,m) , T ,(n,m) / 9 d \ . . 

F 2 V = J 7 ^ - *i - K 2^- 0n, (3.18) 
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^' m) = \i?x - Kq)^t m) + ^4^, (3-19) 



~ 3q / d d \ k\ + k\ d 2 ( d d 

(3.20) 

In (13.151) the Laplace operator is taken only with respect to £, and this problem 
involves s as a parameter. So, we can consider (13. 15j) as a problem for the Dirichlet 
Laplacian S in u with a dependence on s. Since A n is a simple eigenvalue of S, 
the solvability condition of ( 13 . 1 5[) is the orthogonality of the right hand side in the 
equation to (f> n in L 2 (u), 

X (n,m) ¥ n, m) + ( F ^\ ^ ^ = , s G (0,, ), (3.21) 

where we have taken into account the normalization of <p n and the formula (13.81) . 
Let us evaluate the second term in the left hand side of (13.211) . 
Integrating by parts, we obtain 

"" q - K2 £) f-'* = \ I q { k w 2 - K2 4) *- df 

i ( ,2 ( s d \ 4 + 4 

W^ = ^ J K<f> 2 n dt = 0, 
J (PnTZ 2 ^^ = - J \TZct> n \ 2 di = -C n (u), 



k*C b ( W )) ' m) + A^^ H =0, a 6 (0, So ). 



We substitute the identities obtained, (13. 8j) . (13. 18j) . ( I3.20|) into (I3.21[) and arrive at 
the equation 

a^r^ (4+4 

ds 2 V 4 

Together with the boundary condition (13.91) it can be rewritten as 

C n ^' m) = X^' m) ^' m) . (3.23) 

Thus, Aq is a (simple) eigenvalue of C n and ^q 1 '™^ is the associated eigenf unction. 
In what follows we assume the eigenfunctions ^/^'^ are assumed to be orthonor- 
malized in 1/2(0, so). We also note that by the smoothness improving theorems 

^ n ' m) eC°°[0,s o ]. 
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Let V n be the orthogonal complement to {0„} in L 2 (cu). By we denote 
the restriction of S on V n D W / 2 2 (w). It is clear that the operator (<S^ — A)" 1 is 
well-defined in and bounded as that from V n in W$ (uj). 

It follows from (E22D that 



The identity fl3T2Tj) is satisfied due to f l^23|) and it yields X^' m) ^ n,m) +F 2 (n ' m) G 14- 
Hence, by (13.181) . (I3.24p we have (J 7 + \ l ^' m ^)ip^' m ^ G V^,. Taking into account this 
fact, (13.181) . and (13. 12p . we return back to the problem (I3.15p . and write its solution 

as 

# m) (s,0 = V?" m) (s,0 + \^r\s)MO^o + ^ m \s)MO, (3.25) 

where the function \]/ 2 n ' m ) is unknown and should satisfy the boundary conditions 

Bearing in mind the belongings n G C°°(uJ), G C°°[0, So], and the identities 
(13. 20 p . (II. 3p . it is easy to see that the function (J 7 + X^ 1 ' 171 ^)^^'^ can be represented 
as a finite sum 

(? + Ai"' m) )4 n ' m) = E^T^T^), 

where # 2 n j m) G C°°[0, s ], F 2 ( "' m) G C°°(u;) n V n . Thus, 

By the smoothness improving theorems (p^f^ ^ C°°(Zv)! and therefore ip^ £ 

c°°(n). 

We substitute the formulas (JHJ|, dS2D, M, fl3~T0p . (EH, (jXH)) . (TXT7j) . 
fl3T20p . ( 13~25P into the problem ([33]) for z = 3, 



_A f - A n )4 n,m) = \p m) ^' m) + Ai n ' m VS n ' m) + F^ m) in 

4 n>m) = o on an, 

d d 



F 3 (n ' m) :=(Jk - A n q)$ n ' m) + ^ 2 q4 n ' m) + (^3 - A^ m) q)#' m) . 



(3.26) 



(3.27) 



We again treat this problem as that for S and depending on s, and the corresponding 
solvability condition is 

, (n,m) / , (n,m) , \ \ (n,m) T f (n,m) . / j-i(n,m) , \ _ n oo\ 

^0 Iri ^ 9n)L 2 {u) + A x W + (^ 3 , <pn)L 2 (w) ~ U. (3.28) 
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In the same way how the equation (13.231) was derived we obtain 

(Cn - Aj"' m) )^S n * m) = Ai n ' m) ^ n ' m) + ft' m \ (3.29) 

q n (s) := (q(s, -)0n,0n)L 2 H = ^l(s)(60n,0n)L 2 H ~ K 2 (s) (£ 3 0„, 0„)l, 2 ( w ) . 

Since Ag n ' m ' 1 is an eigenvalue of £„, the solvability condition of the last equation is 
the orthogonality in L 2 (0, s ) of its right hand side to the eigenf unctions associated 
with Aq™'" 1 ^, i.e., it should be orthogonal to \& < £ l ' m \ It gives the formula for \^ L,m \ 

A l - -1/3 ' W Ji2(0,s ) 

_ (i?(n,m) ,(n,m)\ ^(«>™)/j,M ™/i( n > m h (3.30) 

--(-^3 jL 2 (fi) - 2 ^° ' q ^° ^2(n)- 

Let us calculate the right hand side of this identity. Integrating by parts and 
employing (EUO]) . (1312)1 . (jOBjl . f l3TT9|) we have 



i((A, + A n )^W^ 



-gU"?^^ >4Vo ; ia(n) (3.31) 
i((^ + A^ ) )4 n ' m) ,q#' m) ) L2(0 ) 



|(Jk - qA n )q^' m) + (^ 2 + A? ' m) )tf' m) , q4 n ' m) 

L 2 (H) 

Employing the eigenvalue equation for n , by direct calculations we check 

(n,m) _ o f 9 8 \ ,(n,m) / 2 , 9,\ n ,(n,m) 



(J 7 ! - qX n )q% ' = -3q [^Iq^ ~ K2 ^J ' K + ^2)^0 
We again integrate by parts, 

(IT \ \ 1 (n,m) <(n,m)\ o I 2 ( ® ® 1 / (n,m) , (n,m) \ 

((J r i-qA n )q^ ; ,q^ )l 2 (^) = ~ 3 ^ ~ K2 ^J ^ '™ J £ 

-((^ + ^)4- m) ,q^ m) ) i2(n) 



=2((K 2 + ^)^ m) ,q^ m) ) 



1,2 (Q) 

L 2 (n) 
L 2 (n)- 



Substituting the identities obtained into (13.31 j) . we arrive at 
((Jk - qA„)^' m) , 4"' m) ) L2( n) = 5 ((^ + 4 + 4 + X^ m) )4 n ' m \ q^' m) ) 



L 2 (0) 

(3.32) 
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It follows from (13.71) that 



We substitute this identity and f )3.32p into f)3.30p and integrate by parts, 



■ } {n,m) i(n,m) 
L 2 (Q) V - ' / L 2 (fi) 

\ J L 2 (n) V / L 2 (n) 

fd q d4 n ^ {n , m) \ ( d^ m \ 



i 2 (H) V / L 2 (f2) 

dq. 



I ^T") / (n,m) i (n,m) \ ( 2t"> / (n,m) , (n^m)^-, \ 



In view of (13.61) . (I3.23P we have 



+ 4)4 n ' m \ ^ m) ) Lm - 2(^S n ' m) , qtf ' m) ) L2{ 



=2((Ajr J - 





^0 I (n,m) 



sr-, I (n,m) i (n,m) \ / 2t>2 i (n,m) , (n,m)\ 
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We integrate by parts employing ( 13.81) . 

o ( n~> i (n,m) \ n ( 9 (n,m) , (n,m) 

-2 I Kk 3 ^ , q-0 o I -2\—K 3 n^ y ',q% 



L 2 (Q) 



9^0 / (n,m) \ / , n,m) , (n,m) 

k 3 Kq, % I - I K 3 —Kip K 

/ 1,2 («) S 
o , (n,m) \ / r\ 

^rO n / (n.m) \ . o I n / fowO 



- /0 t> I (n,m) 1 i o I n / frvO <^ , (n,m) 

I K 3 ^ , ftq^o I +2 1 « 3 ^o > ^q^o 

/ L 2 (n) 
/ L 2 (n) ^ 

900 ^ / (n,m)^-, J / 9q / fn,m) , (n,m) 



«9s ' ru n / V 9s 

1,2 (n) 



i, 2 (n) 



o , (n,m) 



ds 



4^?4(*r>) 2 ) Mn) 

In the same fashion we obtain 

dqdipl n ' m) (n,m.)\ 1 / (n , m ) <9 2 q (n , jm) \ 

0o I =o(^o >7^o • (3.36) 

/ L 2 (fl) 



ds ds ' ru / 2V ru '9s 2 

L 2 (Q) 



Employing the identity lZ 2 q = — q, we have 

f 2n~> I (n,m) , (n,m)^-. \ n / 2n2 i (n.m) , (n,m)\ 

-\k z TI% \% 'TZq) -2[K z n% \qip' ; ) 

\ / Lo(il) 



L 2 (Q) 



2(4n^ m \qn^ m) ) Lm + (4n4 n ' m \4 n ' m) Kq) L2{n) 

n 



1 

1 



)L 2 (n)- 



We substitute the identities obtained and ( I3.34p . ( I3.35p . (I3.36|) into (I3.33j) and arrive 
at the formula (II. 6p for \^( l,m \ 

Let V n , m be the orthogonal complement to {9q } in L 2 (0, so), and £^" m be the 
restriction of £ n to V r „ )m nH / 2 2 (0, s )- The operator (C^ m — \^' m ^ l )~ 1 is well-defined 
and bounded as that from V n ^ m into V n>m D WfoCO; s o)- 

The orthogonality condition (13.301) means that the right hand side of (13.291) is 
orthogonal to \&(™' m \ Thus, / \( ri ' m )\[/| ) ri ' m ) _|_ f^ n ' m ^ e V^ )m . In view of this fact we 
can choose a solution to (I3.29P as 

By the smoothness improving theorems ¥™' m) e C°°[0,s o ]. 

The equations (I3.28p . (13.291) being satisfied, the right hand side of the equation 
in (I3.26P is orthogonal to <p n , i.e., 

By (I3T24D it yields 

In view of this fact we can choose a solution to (I3.26P as 

where the function vj/g 71 '" 1 -* is unknown and should satisfy the boundary conditions 

^ m) (0) = ^ m) (s ) = 0, 

and ip^ l ' m ^ 1 6 C°°(Q). This smoothness is proved in the same way as for i/j^ 1 ' 171 ^. 

The remaining problems for ipi, % ^ 4, are solved in the same way as for i = 1,2,3. 
Namely, the solvability condition of these problems is the orthogonality of the right 
hand side to <p n in L 2 (oj) for each s G (0, s ). In its turn, these condition imply the 
problems for The solvability conditions of these problems imply the formulas 
for Aj. The result of this recurrent procedure is formulated in 

Lemma 3.1. There exist solutions to the problems 113. 5\) given by the identities 
4 n ' m \s,0 = ^ m ^,O + ^& m) (^)q(^O0n(O + ^ (n ' m) (^)0n(O, i > 0. (3.37) 
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The functions ^' m) G C°°(tt), #S n ' m) G C°°[0, s ] read as follows, 
F, (n ' m) = 0, i < 2, 



(3.38) 



i=3 

^( Bim) = _ A (n, m ) r ! ^(n f ) + £ A (n, m )^(n >r ) j ? . ^ ^ 

1 i— 3 

i=o 

where A^" 1 " 1 _\(^ m ) are determined by \3. 8\) . A3.13\) . and 

1 i_1 

,(n,m) _ /j-,(n,m) , (n,m)\ 1 V"~"> , (n,m) / T (n,m) , , (n,m)\ . ^ 1 

A i - ,V>0 )L 2 (Q) - g Z^ A i l*t-i-10rn 0.^0 Jl 2 (Q)' Z ^ L 

i=o 

(3.39) 

Remark 3.1. The formula f!3.38j) involves the term It is zero for z ^ 3, and 

this is why it is absent in (13.18p . (I3.27p . At the same time, it is easy to see that it 
comes from the term (J^ — A^^)^-™'™ in (13. 6ft . when ipi-'™^ is taken in accordance 

with Ezp . 

Given any iV ^ 1, we denote 

AT-2 N 

X^f :=^ 2 An+ EA^' ^ m) («,0 := E £ Vf' m) (^0- (3-40) 

i=0 i=0 

The next lemma follows directly from Lemma 13.11 

Lemma 3.2. The function G C°°(Q) and the number X^^ satisfy the equa- 

tion 

{U E -X^ V£ )^ = ht;™\ (3.41) 
where the right-hand side obeys an inequality 



||,(ra,m)|i ^ rt(n,rn) N-l (1 AO\ 

\\ n e,N \\C k (Ti) ^ U N,k £ > \ AAl ) 

(n,m) 



~e,N \\C k (Q) ^ ^N,k 

with constants C^" ,) > independent of e. 
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We proceed to the justification of the formal asymptotics (13. ip . (13. 2p . We use 
the standard approach based on Lemmas 12, 13 from [T5|. More precisely, we use 
these lemmas in the formulation presented in Lemma 1.1 in P31 Ch. HI, Sec. 1.1]. 
For the convenience of the reader below we give the mentioned lemma. 

Lemma 3.3. Let A : H — >■ H be a continuous linear compact self-adjoint operator 
in a Hilbert space H . Suppose that there exist a real fi > and a vector u G H , 
such that \\u\\h = 1 and 

\\Au — [au\\h ck, oc — const > 0. 
Then there exists an eigenvalue fa of operator A such that 

\fa — n\ < a. 

Moreover, for any d > a there exists a vector u such that 

\\u — u\\h ^ 2ad~ 1 , \\u\\h = 1, 

and u is a linear combination of eigenvectors of operator A corresponding to eigen- 
values of A from the segment \p — d, \i + d] . 

Let e be small enough. Denote ipj?^ := Pe^Velv ■ We rewrite (13.411) as 

p^^p-V^) = X^f^f + p-V^M. (3.43) 

The identities (12.31) . (12.41) . ( 12. 5 p imply that the operator W, £ is self-adjoint and non- 
negative. The same is obviously true for pJ 1 ^ 2 'H £ p7 1 ^ 2 . Let 5 > be a positive 
number. Then the operator A E := (1 + ^p^^'HeP^ 2 ) -1 is well-defined as an 
operator in L 2 (fl), is bounded and self-adjoint, and satisfies the estimate 

\\A £ \\ ^ 1. (3.44) 

As it follows from (12 .4p . for any / G -£/ 2 (^) the function v = A £ f is a generalized 
solution to the boundary value problem 

- 6 ^ 1/2 i2w^w p7l/2v+v = f in ^ v = on m - 

i,j=l 1 3 

Hence, the operator A E is also bounded as that from L 2 (fi) into VK 2 x (fi). In view of 
the compact embedding of Wl(Vt) in L 2 (f2) the operator A £ is compact as that in 
L 2 (fi). 

We rewrite (13.431) as 

(1 + s^r 1 ^ = A £ ^ N m) + sh<$?\ h^f : = (1 + ^x^-u^- 1 /^^). 
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It follows from (I3.40p . (I3.42p , (I3.44p that we can choose 5 = 5(s,N,n,m) so that 

I i~(n,m) II 

1 < 1 + 5\ft ] ^ 2, \f " £a(n) ^ C^ m) 5e N - 1 ^ -, (3.45) 

£ i n ' ii , (n,m) ii JV o' 



where Cjy '"^ are some constants independent of e. 
We apply Lemma [3.31 with 

H = L 2 (Q), A = A £ , /I = (1 + SX^' 1 



V,7V 

u = — —, — \ , a = C\r 'e o, 

<• I (n,m) ii JV 

% N \\L 2 (n) 



and conclude that there exists an eigenvalue /i^^ °f A £ such that 

l/4X } " ( X + ^TV 1 ! < Cp m) 5e N -\ (3.47) 

where C^ ,m ^ are some positive constants independent of e. It is clear that := 
- l)^" 1 is an eigenvalue of p e ~ 1/2 ?4p7 1/2 . It follows from (1X15]) . flX47|) 

that 

^ l -T—^ ~ Cp m) e N ' l 5 1 + 5\ { ;f > 2. 

The last inequality and (13 .45 p . (I3.47P imply 

|(1 + SXffi) - (1 + 5A^ m) )| < Cp m) 6e N ^\l + 5A^ m) l |1 + 

|A^ m) ) " ^> m) | < AC^ m) 5e N ~ 1 . (3.48) 

By e^' m ' ) we denote a monotonically decreasing (in N) sequence such that 



r (n,m, < r (n,m, < < 



Letting 



A^ m ) : =A^ m) for [ep m \e (n ' m) 



N i °N+l 7' 

( n > m ) n f ^-V2ay -1/2 



and employing (I3.48p . we see that the eigenvalue Ae ,m of p e H e p e has the 
asymptotic expansion (11.51) . To complete the proof it remains to note that the 
eigenvalues of p £ 1 ^ 2r H £ p £ coincide with those of H. £ . 



4 Proof of Theorem 1.2 



We begin the proof with the result of Theorem 4.4 in [3]. Namely, the item (ii) 
of this theorem says that given any M > 0, there exists Eq(M) > so that for 
all e < e the first M eigenvalues A m (e), m = 1,...,M, of H £ taken counting 
multiplicities satisfy the asymptotics 



A m (£)=^A 1 + A^ m) +o(l) 



m 



M. 



(4.1) 
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Since the eigenvalues Ag 1,m ' ) of C\ are simple, the same is true for the eigenvalues 
A,n(e). 

It follows from (14. ip that there exists a fixed number 9 > so that for each 
m = 1, . . . , M, and all e < Sq{M) the interval 

(e~ 2 X 1 + A? ,m) - 0, £- 2 A! + A^ m) + 0) (4.2) 

contains exactly one eigenvalue of % E which is X m (e). In accordance with Theo- 
rem [ITT] the eigenvalues Ai 1,m ^ of "H e satisfy the same asymptotics as X m (e). Hence, 
for sufficiently small Eq(M) and all m = 1, . . . , M each of the intervals (14. 2 p contains 
the eigenvalue A^ 1,m ^(e). Therefore, A( 1,m )(e) = X m (e), and it proves the statement 
of the theorem on the eigenvalues. 

To prove the statement on the eigenfunctions, we adopt the same notations 
as in the proof of Theorem 11.11 We again apply Lemma 13.31 with (I3.46P and we 
take d = e N ^ 2 . Then there exists a linear combination ipf'^ of the eigenfunctions 
associated with the eigenvalues of A e lying in the segment 

[(1 + SX^fy 1 - e N/ \ (1 + SX^fy 1 + e N/2 ] (4.3) 

such that 

|L7(l,"i)||„l/2„ / ,(l,m) II -1/2 ,/,(!."») II ^ Mm) N/2-1 r (a a\ 

\\%,N llPe Ve,N l|l*(n) ~ P e ' ^ e ,7V || L2 (fi) ^ G N £ °' I 4 - 4 ) 

where are positive constants independent of e and 5. We choose 5 = e 2 and 
by the asymptotics (11.51) for Xi 1 '^ we obtain the estimate 



1 1 



1 + 5X?' m) 1 + 5X ( '' P) 



sixi 1 '^ - xi 1,p) \ 



,(1,™)| | "I _l §\0->P) 



>Ce 2 



for e small enough, m,p = 1, . . . , M, m ^ p, where C is a positive constant in- 
dependent of e, m, and p. Hence, for N ^ 5, m = 1, . . . , M, and £ small enough 
the intervals (14 .3p contain no eigenvalues of A £ except (1 + e 2 \^' m ')~ l . This eigen- 
value is simple, since the corresponding eigenvalue X^'™^ is simple. Thus, the linear 
combination ipf'^f 1 is an orthonormalized in L 2 {VL) eigenfunction associated with 

(1 + £ 2 Ae Moreover, it is independent of N. 

By the defintion of ipe ' and Lemma 13.11 we have 

N 
j=0 

for all N ^ 0, m — 1, . . . , M. Hence, there exists a function c m = c m (e) such that 

c m {e) = \\ V \l 2 ^f\\ L2m + 0{e N+1 ) for all N > 0. 
The identity obtained and (14.41) yield 

mr ] - v'MT ikw = 0(^ /2+1 ), := Cm (£)^ m) - 
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Denoting ?/4 (s,£) := p £ (s, £)i/Je,* m (s, £), we can rewrite the last equation as 
ll^lk(n)=0(^ 1 ), ^:=Vf' m) -^- (4-5) 
It is also clear that the function ipi 1 solves (12. 6p with A e = Ae ' . Hence, by 

= A^p^ + hW, h%> := (A&"> - A^) Pe ^ - h?f. (4.6) 
Due to this equation we can write the integral identity 

{H £ <P £N , <S> £tN )L 2 {n) - K >{Pe® £ ,Ni ^e,N)L 2 {n) + {n, £jN , ® £:N )L 2 (n)- 
From (Oj) . (JS32D, and (JTSJ) we derive 

ft#llo<n) = 0(e^ +1 ), iV^5, fc^O. (4.7) 
Together with ([231), (D, (@~5D it gives 

II V W )^ll! 2( n) < CCAWV^^, V (s>0 ^) L2( n) = 0(^" 4 ), m — 1, . . . ,M. 

(4.8) 

Combining (I4.5p and (14. 8p . we conclude that the asymptotics (11.71) hold true in 
IV^ (fi)-norm. 

We proceed to the proof of (11.71) in C fc (f^)-norm. First we note that by the 
standard smoothness improving theorems we have ipi 1,m ^ G C°°(fK*)) for all t G 
(0, s /2). The rest of the proof follows from flOji . (H77j) . (Oj) . the embedding 

of W 2 fc+2 (fJ) into C fc (fi), and from the next lemma applied to 

Lemma 4.1. Let u s be a solution to the equation 

'H £ u £ = X £ 1 ^p £ u £ + h, heC°°(T}). 



Then u £ G C°°(fiM) for all t G (0, s /2), and 

IHL*(nW) ^ Cfe£ _2(fe-1) (lk £ ||w 2 i(n) + ll^ll w *- a (n)) ( 4 - 9 ) 

for all k ^2, t G (0, Sq/2), where the constants are independent of e, h and u £ . 

Proof. The smoothness of u £ follows from that of , h, p e , and the smoothness 
improving theorems. For the sake of brevity we denote £1 := s, £ := (£1,^2, £3)- 

Let = Xi\^i) be an infinitely differentiable cut-off function equalling one 
as £1 G [t, s - t] and vanishing for & G [0, t/2] U [s - t/2, s ]. We fix i G (0, s /2) 
and denote 

It is straightforward to check that Ue solves the equation 

U £ uf = Afc-Oprf + /i X W + L £, u £ , |g) . (4.10) 
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o ft) 

The symbol ^ in the arguments of G\ indicates that this function depends on all 

first derivatives of u £ . The function G® is linear with respect to u £ and the first 
derivatives of u £ . The coefficients at u £ and at the first derivatives of u £ belong 
to C°°{p) and are bounded in C K (f2)-norms uniformly in e for all k ^ 0. We 
differentiate (I4.10p with respect to £1, 



du® 



P- 



dv^_ 







m 1 hx; 



du 



s~ 2 G® (e^,u £ 



d 2 u® 
dm. 



(4.11) 



The symbols q^ U q^. i n the arguments of G 1 ^' indicates that this function de- 

pends on all first and second derivatives of u® . The function G® is linear with 
respect to u £ and the first derivatives of u £ . The coefficients at u £ and at the first 
derivatives of u £ belong to C°°(JY) and are bounded in C K (f2)-norms uniformly in e 
for all k ^ 0. 

Starting from (14. lip and proceeding as in (14. 8p . we obtain 



At) 




® „du®' 
V- 



L 2 (Q) 



(0 



e' 2 G\ 



(t) 



,(*) 



96 ^i^- 
1 



L 2 (n) + ll w e|lwi(n) J + 



L 2 (n) 



d 2 uf 
du® 



V 



du® 



L 2 (Q) 



2 

L 2 (Q)' 



(4.12) 



Here and till the end of the proof by C we indicate non-specific constants indepen- 
dent of e, u £ and h. The obtained estimate implies 



.du® 



L 2 (n) 



£ Ce 



\h\ 



L 2 {Ci) 



\u 



(4.13) 



Let us estimate 



' du. 



(t) 



L 2 (0) 



2,3. We could have tried to differentiate 



(14. lUp with respect to £j and proceed as above. However, the function does 
not vanish on (0, sq) x du and this is the main difficulty This is why we have to 
employ a slightly different trick. We introduce an infinitely differentiable in u cut- 
off function \2 = X2(0 equalling one in a small neighborhood of du and vanishing 
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outside a bigger neighborhood. Writing the equation for Ue (1 — X2) similar to 
f)4.10p and proceeding as in (14.111) . (I4.12|) . one can show that 



(4.14) 



In a small neighborhood of duo we introduce new variables ( = (0, (2, C3)) where 
Ci = £1, C2 is the arc length of du and (3 is the distance from the point to duj 
measured in the direction of the inward normal. Then it follows from (I4.10p that 
the function = X2^P satisfies the boundary value problem 



f)2jt) 3 ~ (t) 

6 V * _/ )v v (*) j_ F -2 V p (t) R (e)..(t) , V- R {e)OUe 

CG(0, So )x(0,4 0) )x(0,Cf), 



i=l 



(4.15) 



where is the length of £g U; is a small fixed number. The operator in the 
left hand side of the last equation is elliptic uniformly in C and e. The function v, 



.(0) 



(t) 



satisfies periodic boundary condition as (2 = and (2 = C2 > anc ^ vanishes as £3 = 0, 
C3 = (3°^ Ci — 0, Ci = s o- The coefficients B^\ are infinitely differentiable and 
satisfy the estimates 

\\B<f? 



ij »Cl(|P,«o]x[o4 0) ]x[0 J Cn) 



>(0), 



LB; 



00 I 



Ic^acMolxio^'lxio.cl ']) 



^ Ce 



We differentiate (14.151) with respect to £2, and in the same fashion as in (14.111) . 
1Q2]) . (jUSD we obtain 



^2 



L 2 ([0,so]x[0,c( 0) ]x[0,cl 0) ]) 



^ Ce 



(4.16) 



Now we express the term 
gives the estimate 



from ( jUSJ) . Together with (Q5j» . gUD, P~TB]) it 



<9Cf 



L 2 ([0,s ]x[0,^ 0) ]x[0,ci 0) ]) 



^ Ce 



\h\ 



L 2 (fi) + ||W e || W2 i(f2) ) • 



This estimate and (14.131) . (I4.14p . (I4.16P imply (14. 9p for k — 1. To prove it for other 
fc's, it is sufficient to proceed as above starting with differentiating the equations 

2,3. □ 



iUI SCI ' dQ^ 1 



\ (t) s (*) • 
X2K , ^"X2«e , « 
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